Hydrodynamic effects in driven soft matter 
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Recent theoretical works exploring the hydrodynamics of soft material in non-equilibrium situ- 
ations are reviewed. We discuss the role of hydrodynamic interactions for three different systems: 
i) the deformation and orientation of sedimenting semiflexible polymers, ii) the propulsion and 
force-rectification with a nano-machine realized by a rotating elastic rod, and Hi) the deformation 
of a brush made of grafted semifiexible polymers in shear flows. In all these examples deformable 
polymers are subject to various hydrodynamic flows and hydrodynamic interactions. Perfect stiff 
nano-cylinders are known to show no orientational effects as they sediment through a viscous ffuid, 
but it is the coupling between elasticity and hydrodynamic torques that leads to an orientation 
perpendicular to the direction of sedimentation. Likewise, a rotating stiff rod does not lead to a net 
propulsion in the Stokes limit, but if bending is allowed an effective thrust develops whose strength 
and direction is independent of the sense of rotation and thus acts as a rectiflcation device. Lastly, 
surface-anchored polymers are deformed by shear flows, which modifles the effective hydrodynamic 
boundary condition in a non-linear fashion. All these results are obtained with hydrodynamic 
Brownian dynamics simulation techniques, as appropriate for dilute systems. Scaling analyses are 
presented when possible. The common theme is the interaction between elasticity of soft matter 
and hydrodynamics, which can lead to qualitatively new effects. 



I. INTRODUCTION 

Dynamics in soft matter systems far from equilibrium 
is of central interest in various interdisciplinary fields 
like biophysics or microfluidics. Most generally, soft 
objects subject to thermal noise show thermally ac- 
tivated shape fluctuations and are sensitive to even 
minute external forces. Moreover, the dynamics of 
small objects evolving in liquid solvents, e.g. water, 
is modified by the presence of flow and viscosity. As 
an example, external forces exerted on micro-objects 
lead to long-ranged effects due to momentum diffusion 
and modify any motion and flow. In particular, dif- 
ferent moving objects or different moving parts of a 
macromolecules are coupled among each other via the 
hydrodynamic flow fields emanating from the differ- 
ent regions. Such hydrodynamic interactions, mediated 
by the viscous solvent, are well studied in the context 
of colloids [IHl] and polymers [S]. They are respon- 
sible for many striking features such as propulsion at 
the micron scale or measurable deviations from pure 
Brownian random dynamics in sedimentation and elec- 
trophoresis experiments [5]-[7|- For small objects and 
rather modest fiow velocities, the non-linear term in 
the Navier-Stokes equation is negligible and one ob- 
tains the Stokes equation which describes the so-called 
creeping flows. Since the hydrodynamic equations gov- 
erning the flow become then linear, so are all response 
relations between forces and torques applied on rigid 
micro-objects and the resulting velocities [I]. This lin- 
ear relationship is completely defined by the shape of 
the objects and thus illustrates the importance of shape 
design in microfluidics, for instance. However, the is- 
sue becomes more complicated for deformable soft ob- 
jects like membranes, soft colloids or elastic rods and 
macromolecules (cylindrical viruses, synthetic or bio- 
polymers). In that case, the deformation of the object 



depends on the flow which thus changes hydrodynamic 
boundary conditions in an intricate and non-linear fash- 
ion, since the object shape acts back on the flow field via 
hydrodynamic screening. Solving the problem becomes 
much more complicated, both at the numerical and an- 
alytical level. Some analytical results have been ob- 
tained either close to equilibrium where linear response 
theory applies or by invoking the pre-averaging approx- 
imation, as used in the context of the Zimm-model for 
dynamics of polymers '5\. However, to unravel the full 
coupling between shape and flow in situations far from 
equilibrium, simulation techniques prove very useful. 
In this review, we focus on three situations involv- 
ing hydrodynamics of soft matter far from equilibrium: 
i) the sedimentation and orientation of elastic rods, 
ii) the propulsion and force-rectiflcation by a rotat- 
ing semi-flexible filament, and Hi) the deformation of 
a brush of grafted semiflexible polymers in shear flows. 
As will be discussed, the selection of these topics is 
closely related to novel emerging experimental tech- 
niques and activities. 

Example i) is related to the effect of anomalous bire- 
fringence, which is observed in experiments involving 
charged rod-like particles such as Tobacco-Mosaic [5], 
FD viruses |9j or synthetic polyelectrolytes [101 E] and 
cylindrical micelles [T^ in electric fields. Normally, 
charged objects in electric fields orient parallel to an 
externally applied electric field due to the anisotropic 
polarizability which is mostly caused by the easily dis- 
placeable counterion cloud. This is called normal bire- 
fringence. For almost all rod-like charged systems, one 
finds in a certain parameter range also anomalous bire- 
fringence, where particles orient perpendicular to the 
direction of the electric field and motion. This effect 
is obtained for long particles, law salt concentration or 
particle concentrations beyond mutual overlap. A sim- 
ilar effect, namely perpendicular orientation of moving 



cylindrical objects, is in principle also observable with 
sedimenting rod-like particles [13j . This phenomenon 
is, at present, not completely understood [Hj. A pos- 
sible relation with a mechanism for hydrodynamic ori- 
entation associated with rod deformation was recently 
suggested [15]. As follows from the linear-Stokes equa- 
tion, perfectly straight and stiff nano-cylinders show no 
orientational effects as they sediment through a viscous 
fluid. But an elastic rod is deformed due to hydrody- 
namic forces, and the coupling between elasticity and 
hydrodynamic torques leads strikingly to an orientation 
perpendicular to the direction of sedimentation. 

In example ii) the question of propulsion of a microm- 
eter sized object through a viscous fluid is addressed 
from an engineering point of view. Due to the linear- 
ity of hydrodynamic equations, propulsion on the nano- 
meter scale calls for design strategies very different from 
the macroscopic world, since inertia plays no role and 
friction is the only way of producing thrust |16| I17j . 
The studied propulsion mechanism is inspired by bi- 
ological systems such as bacteria, which move by ro- 
tating propeller-like appendices, so-called flagella [TSL 
[20] . Those propellers involve helical stiff polymers that 
are rotated at their base by a rotary molecular motor. 
Hydrodynamic friction converts the rotational motion 
of the helix into thrust along the helix axis 16, 18J. 
For a number of biomedical applications, e.g. for di- 
rected motion of artificial viruses through cells or nano- 
devices through the bloodstream, it is desirable to de- 
velop similar synthetic propulsion mechanisms or to in- 
corporate biological single-molecule motors into syn- 
thetic environments. A second possible field of ap- 
plication involves mixing strategies in nano-fluidic de- 
vices, which could be achieved by moving or rotat- 
ing surface-anchored polymers [21] . Recent discoveries 
opened the route to the synthetic manufacture of ro- 
tary single-molecule motors driven by chemical [22^ or 
optical |23l [24] energy. An ATPase can also be fixed to 
different substrates and used to rotate metallic [25] or 
organic nano-rods [35]. All these experiments raise the 
question about the minimal design necessary to con- 
vert the rotational power of such nano-engines into di- 
rected thrust in a viscous environment. Since stiff he- 
lices are actually difhcult to manufacture down to the 
micron scale, we investigate whether a straight elastic 
filament can act as a propeller. Our design consists of 
a straight elastic rod that rotates around a point with 
a constrained azimuthal angle, i.e., on the surface of 
a cone. In the absence of elasticity, no effective thrust 
is produced as one averages the motion over one whole 
cycle of rotation. But a rod with finite bending mod- 
ulus will be deformed and in fact take on a shape that 
resembles one period of a helix. Due to the symmetry 
breaking, a net thrust is obtained which in fact does not 
depend on the sense of rotation of the rod: the device 
acts as a nano force- rectification device [57] . 

Finally, we study the infiucnce of a shear flow on the 
conformations of surface-anchored semiflexible poly- 
mers. When the shear rate increases, polymers will 
be oriented in the direction of the shear and be bent 
down in order to be more and more parallel to the sur- 



face. These conformational effects clearly depend also 
on the grafting density of polymers and on the stiffness 
of polymers. As polymers bend down onto the surface, 
the hydrodynamic boundary condition changes and the 
flow can penetrate closer towards the surface. Again, it 
is the coupling between hydrodynamic fiow effects and 
elastic deformation of polymers which leads to intricate 
non-linear effects. These phenomena have been exper- 
imentally investigated using colloidal spheres with sur- 
face anchored DNA molecules, which were being held 
fixed using a laser trap in uniform flows of different 
strengths 28 . Since the force acting on the sphere 
is measured in such an experiment, the hydrodynamic 
drag and therefore the position of the surface of shear 
can be measured with high precision. Similar phe- 
nomena also appear in the glycocalix layer, which is a 
polymer-brush like coating of the endothelial cell layers 
in blood vessels. The shear rates in blood stream can 
be enormous, and conformational changes of the glyco- 
calix layer has important consequences on permeability 
of drugs, nutritional agents or viruses. 
The paper is divided in five sections. In Section [iT] we 
give an introduction of hydrodynamics at low Reynolds 
numbers and focus particularly on hydrodynamic inter- 
actions between two moving spheres in a viscous fluid. 
The integral representation of the flow velocity, intro- 
duced by Oseen using Green's function techniques, is 
developed, since it is at the heart of our simulation 
methods. We present Green's tensors both for an un- 
bounded fluid and a semi-infinite fluid medium with a 
no-slip wall. Section [fll[ presents our model of semi- 
flexible polymers under external force which is imple- 
mented in Brownian dynamics simulations in the pres- 
ence of hydrodynamic interactions. Section |V| is de- 
voted to hydrodynamics of elastic polymers under ex- 
ternal force load. After recalling some basic results for 
rigid cylinders and helices (Section IV ) , we focus on the 
coupling between hydrodynamic interactions and elas- 
tic deformations. This coupling leads to the orientation 
of elastic cylinders under an homogeneous field, propul- 
sion with rotating polymers, and deformation of brush 
layers in shear fiows. 



II. HYDRODYNAMICS ON THE NANOSCALE 

A. Hydrodynamics at low Reynolds number 

The Navier-Stokes equation combines both convective 
and viscous terms, the ratio of which is given by the 
Reynolds number. Re = - — , where U and L are the 
typical velocity and size of the flow, p the mass per 
unit volume and rj the shear viscosity of the fluid [29] . 
In this paper, we are dealing with small objects mov- 
ing slowly in a fluid, like a polymer, a macromolecule 
or a latex bead. The typical scales are in the range 
L ~ lO-^-lQ-^m and U ~ 10-**~10-''m/s yielding in 
water (p = lO^kg/m^, 77 = lO^^p^/g) ^g ^ iQ-iO- 
10^^. Hence, hydrodynamics on the nanoscale be- 
long to the small Reynolds number regime. Re <C 1. 
In this framework, viscous damping (momentum diffu- 



sion) dominates inertia (momentuni convection), and 
the creeping flow produced by moving objects is de- 
scribed by the Stokes equations [1, i29i 
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where p{r) is the pressure field, v(r) the fluid velocity 
and f a point force at r = 0. In this paper, the fluid 
can be taken as incompressible at the time scales we 
consider. 

The solution of the linear equations (fl )-([2]) can be found 
using Green function techniques |30H32[ . The Oseen 
tensor i33, , also called Stokeslet, G(r) describing the 
fluid velocity disturbance caused by a point force ex- 
erted at the origin, is 
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where r — |r| and r (g) r = rr"^ is a dyadic tensor. The 
pressure field is the potential of a dipole distribution of 
strength — f . Since we are dealing with an unbounded 
fluid, the flow disturbance vanishes at inflnity: u(r — >■ 
oo) =0. The famous inverse first power dependence 
in r means that hydrodynamic interactions are long- 
ranged. Moreover the flow velocity v(r) is not uniform, 
the perturbation being larger in the direction of f . This 
non-uniformity is linked to the symmetry-breaking into 
longitudinal and transverse modes in fluids. 
The flow around any moving rigid body (volume V 
and surface S) in an unbounded medium can then 
be theoretically constructed by using a distribution of 
Stokeslets on the body surface. Indeed, if there is no 
source of momentum in the fluid, the flow is only de- 
termined by the no-slip condition at the body surface. 
For a steady creeping flow, the propulsive external force 
applied to a moving object is balanced by the frictional 
force exerted by the fluid on the body surface, which is 
given by integrating the stress tensor 



o- = -pi + r/[V (g) V + (V «) v)^] 



(5) 



over the body surface [2U]. Most generally, the bulk 
fluid velocity is then given by the surface integral 



v(r)= f Gir-r')-i{r')d^r' 
Js 



(6) 



where the surface force distribution f (r') is the source 
of the flow fleld. 

As an illustration of these concepts, we calculate in the 
following the flow induced by a moving sphere in an 
unbounded fluid. It will allow us to introduce hydrody- 
namic interactions between spherical particles evolving 
in a creeping flow. 



B. Flo'w around a moving sphere and 
hydrodynamic interactions 

For a rigid object with an unspecifled shape, the calcu- 
lation of the flow fleld at distance r from the immersed 
object of size a using Eq. (|6| becomes hardly tractable 
and, similarly to electrostatics, a multipole expansion 
can be done in powers of a/r. At flrst order, one can 
approximate the velocity distribution by a point force 
exerted at the centre of the body. This corresponds to 
the far fleld solution, which is a Stokeslet. However, 
this approximation is no longer valid when we consider 
the flow fleld closer to the immersed body. We shall 
then take into account its flnite size. Let us consider 
a rigid sphere of radius a with its centre positioned at 
Ti. We rewrite Eq. m using a Taylor expansion and 
obtain the fluid velocity Vj (r) , associated to a uniform 
force density ii{r) — fj applied on the sphere surface 
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where F^ = ^gfi(r')d r' — Ana^ti is the total force 
acting on the particle j. This result is indeed exact for 
a sphere since its high symmetry ensures that higher 
order surface integrals are zero. 

One can now calculate the velocity of the rigid sphere, 
v,f, by utilizing the no-slip boundary condition on the 
sphere surface: the sphere velocity is equal to the fluid 
velocity at the surface. We thus flnd the Stokes law 
for the velocity of a sphere which moves slowly under a 
force F, as [DIM] 
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Since the sphere is moving with a stationary velocity, 
— Fi is also the frictional force acting on a sphere mov- 
ing at velocity v|. 

Suppose now that we have two spheres i and j of same 
radius a in an unbounded fluid with two external forces 
Fi and F^ acting on them. The fluid velocity gener- 
ated by Fj on particle j also tends to drive the sphere 
i: these interactions via momentum diffusion are called 
hydrodynamic interactions. Of course, the modifled ve- 
locity of particle i should in return modify the flow fleld 
close to particle j and then its velocity and so on. How- 
ever, these effects are successively smaller by a factor 
2a/rij where Vij = r^ — Yj and can be calculated by the 
method of reflections [U [35] . In the following, we con- 
sider that the spheres are widely separated, 2a/rij <^ 1, 
and the velocity of particle i is simply the sum of two 
terms: its self-propulsion by F^ given by Eq. (Is]) and 
the contribution of the flow created by the sphere j, 
Vj (r) , given by Eq. (l7| . The sum of the two effects can 
be written as 
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In the above equation we use the Taylor expansion of 
the flow field with respect to the particle centre r^ . This 



naturally leads to the mobility tensor, that describes 
hydrodynaniic interactions between particles, as 
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where 1 is the 3x3 unit matrix. One notes that this 
mobility tensor includes the lowest corrections of par- 
ticle size over the Oseen tensor description, and cor- 
responds to the well-known Rotne-Prager tensor in an 
unbounded fluid [36]. One thus should keep in mind 
that our treatment of hydrodynamic interactions is only 
approximate for particles which are separated by a dis- 
tance on the order of a. 



Hydrodynamic interactions near a no-slip 
surface 



If a particle is moving near a rigid boundary like a solid 
planar surface, the situation becomes much more com- 
plicated due to boundary conditions. The fluid veloc- 
ity must vanish at the surface which is at rest. Here 
again, theoretical methods used in electrostatics apply. 
Blake [37] developed the image method for hydrody- 
namics near a plane and obtained the Green's function 
in the presence of a no-slip planar boundary at z = 
as 

Gw/(r„r,) = G(r)-G(R)+G^(R)-G^^(R) (12) 



where r 



r, and R 



: — Tji and j' is the im- 
age with respect to the plane z = (see Fig. [T]). The 
Blake's Green function relates the flow at r^ to a unit 
point force applied at Vj in the bounded fluid medium. 
One notes that the no-slip boundary condition on the 
plane of the wall is satisfied by the image system [551 - 
|44] . consisting of the original Stokeslet G(r), the im- 
age Stokeslet G(R), a Doublet G^(R) and a source 
Doublet G'^-'^(R). If one uses the tensor G^aii given 
by Eqjl2]), instead of G'^ given by Eq. (13]), in equa- 



tion (11[), we obtain the equivalent of the Rotne-Prager 
mobility tensor in the presence of a no-slip wall, which 
is the level of description we chose for the last example 
of surface-anchored polymers in shear flows. 



III. HYDRODYNAMIC NON-EQUILIBRIUM 
SIMULATIONS OF SOFT MATTER 

A. Review of numerical methods 

On the methodological side, the adequate inclusion 
of hydrodynamics in a theoretical description of soft- 
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FIG. 1: Top: Two spherical particles i and j, where r = 
Ti — Tj, in the presence of no-slip planar surface. The dashed 
sphere j' beneath the surface represents the image of j-th 
particle and R = r^ — r^/. Bottom: Schematic picture for 
determining self-mobility of a single particle where R is the 
relative position vector between the particle and its own 
image. 



matter systems has been a long-standing problem |45l - 
W7\ . From a simulation point of view, molecular dy- 
namics (MD) simulations with explicit consideration 
of solute and solvent take into account hydrodynamic 
interactions on an adequate level. However, most of 
the simulation time is used to evaluate the solvent 
dynamics in great detail, which is not necessary to 
achieve the proper hydrodynamic behavior on the col- 
loidal or polymer length scale: indeed, MD simulations 
require time-resolution on the order of 10~^^s to de- 
scribe inter-molecular forces whereas the typical diffu- 
sion time of a colloid of size a = lOnm in water is 
r = Qurjo? /kBT ~ 5 x 10^^ s. Moreover, it will displace 
approximately 10* solvent molecules. This large length- 
and time-scale gap between the solvent molecules and 
the solute calls for a coarse-grained and simplified de- 
scription of the solvent dynamics. 

This has led to the development of novel mesoscopic 
simulation techniques. Prominent examples are lattice 
models such as lattice-Boltzmann methods [^5H5T] , par- 
ticle based off-lattice methods such as dissipative par- 
ticle dynamics [52H54] . and multi-particle-collision dy- 
namics [SS] [S5] , and more coarse-grained scheme such 
as Lowe- Anderson thermostat [ST] [55] and stochastic 
rotation dynamics [59j . 

We next describe the implementation of a Brownian dy- 
namics simulation for elastic polymers. We first intro- 
duce the general framework of Langevin dynamics with 
hydrodynamic interactions and then adapt it to semi- 
fiexible polymers. In contrast to lattice-based methods, 
for which the numerical complexity scales with the sys- 



tern volume, hydrodynamic simulation methods |601 - 
I62j lead to a numerical complexity which scales with 
the number of monomers included in the simulation. 
Clearly, hydrodynamic simulation methods will be ad- 
vantageous whenever dilute systems are being stud- 
ied and the number of polymers is small. Likewise, 
semi-infinite systems {e.g. polymers grafted or ad- 
sorbed to a single solid surface and connected to bulk 
solvent) can be straightforwardly simulated using hy- 
drodynamic simulation methods and the appropriate 
hydrodynamic Greens function satisfying e.g. the no- 
slip boundary condition at the surface. Following this 
scheme, the solvent is treated as a continuum, which is 
acceptable for flows on length scales larger than the size 
of solvent molecules [17] . Stokesian dynamics does not 
allow to describe the real flow field in the immediate 
vicinity of our structural units, the monomers, whose 
size ranges between 1 — 10 nm. On the other hand, 
the flow field calculated far away form the monomers 



by using the Rotne-Prager tensor, Eq. ( 11 ) remains ex- 
act at the order of (a/r)'^ and is not influenced by the 
flow close to the monomers. Moreover, the hydrody- 
namic interactions between two monomers in contact 
are negligible compared to direct monomer-monomer 
interactions (stretching and bending) or excluded vol- 
ume effects. Hence, our treatment of hydrodynamic in- 
teractions between monomers far apart should remain 
valid. This is borne out by recent Non-Equilibrium 
Molecular Dynamics Simulations of sedimenting carbon 
nanotubes in water, where it was demonstrated that 
the Stokes equation describes simulated solvent flow- 
proflles down to the sub- nanometer length scales |63j . 
More subtle is the question of the appropriate hydro- 
dynamic surface-boundary condition, as it depends on 
surface details and displays variable amounts of slip at 
hydrophobic surfaces [63j . In our case we are mostly 
concerned with hydrophilic objects, where the no-slip 
boundary condition is valid to a very good degree. 



B. Brownian hydrodynamics 

The two following assumptions provide a frame of ref- 
erence for the description of Langevin dynamics 

• the mesoscopic size of polymers and colloids (a '^ 
10 nm to 1 /zm) ensures the validity of Langevin 
dynamics. Fluid particles are then treated as 
variables having fast dynamics and are integrated 
out, whereas material particles follow slow dy- 
namics and are treated explicitly; 

• we assume local thermodynamic equilibrium,, i.e. 
particle velocities obey fast dynamics which im- 
plies that the memory of velocities is lost much 
quicker than the time scales of interest. We are 
thus interested in time scales larger than the mo- 
mentum relaxation time, t ^ T,n — m/{6ari), 
which constitutes the diffusive regime. 

Our polymers are coarse-grained at the nanometer scale 
and modeled as an assembly of M -1-1 spherical particles 
submitted to external and interaction potentials. In 



this framework, we are able to write the time evolution 
of all sphere positions ri{t), which is governed by the 
position-Langevin equation 64n69j 
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+ fcsT^Vr, -/x.^+^.W (13) 

which relates the velocity of sphere i to forces applied 
to all beads including hydrodynamic interactions as 
given by Eq.Q. One notes that the hydrodynamic 
mobility tensor /.tj introduces long-ranged interactions 
and couples distant particles. The stochastic term, the 
Langevin random displacement ^(t), mimics the ac- 
tion of a thermal heat bath and obeys the fluctuation- 
dissipation relation 



(|.(i)®|,(t')>-2fcBT/x„-<5(t-t') 



(14) 



which is numerically implemented by a Cholesky de- 
composition [SSj . We neglect all memory effects in this 
work. In the most part, we study the dynamics of elas- 
tic rods in an unbounded fluid. In this case we use the 
mobility tensor /x^ ■ given by Eq. (11 1 and it is straight- 
forward to check that Vf^. • /Xj = U for every bead and 
the gradient-mobility term in the Langevin equation 
vanishes identically. This is different in the presence of 
a no-slip surface, where this term is non-zero and has 
to be taken into account. 

The forces have two contributions: the internal forces, 
which are derivatives of the interaction potential U, and 
the external force, F°'^', applied on bead j. These forces 
depend on the specific details of the studied system 
and will be described in the following sections. Hydro- 
dynamic interactions between spheres are implemented 
via the Rotne-Prager mobility tensor given by Eq. (fTTl) 
for an unbounded medium and by using the procedure 
defined in Section III CI for a semi-infinite medium with 
a no-slip wall at z = 0. It ensures that the mobility 
matrix is positive-definite which is necessary to use the 
Cholesky decomposition j68j . 



For the numerical iterations, we discretize Eq. ( 13 ) with 
a time step A and rescale all lengths by the sphere 
radius a according to f^ = r/a. The iterative Langevin 
equation in terms of the discrete time variable n = t/A 
now reads 



-VtV{n)+Ff{r 
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M 
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3=0 

where U — U/ksT is the dimensionless potential, 
F^''* — aF^^^/ksT the rescaled external force, and 
the rescaled random displacement has variance unity 
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We express all mobilities 



in dimensionless form by defining /ij — A/Xj /csT/a^. 



The rescaled bare mobility 



Amo 



(16) 



is the diffusion constant in an unbounded space in units 
of the particle radius a and time step A. 



C. Semi- flexible polymer model 

An elastic polymer is modeled as a chain of connected 
spheres (monomers) that interact via stretching and 
bending forces. The non-dimensional elastic potential 
U is the discrete version of the extensible worm-like 
chain model 



U{{rk}) = 
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(17) 

where 0i is the angle between neighbouring bonds of 
sphere i. The parameters 7 = ja/kgT and i — 
e/^aksT) are the adimensionnal stretching and bend- 
ing moduli. For an isotropic elastic cylinder with radius 
a, bending and stretching moduli are determined by 
Young's modulus Ey as e = EyTra'^/'i and 7 = Eyna'^, 
which gives the relation 7 = 4e. If not stated other- 
wise, we always consider isotropic elastic rods in the 
following. The persistence length, i.e. the correlation 
length of tangent vectors along the polymer contour, is 
given by £p = e/{kBT). Twist and torsional degrees are 
omitted since for most synthetic polymers, free rotation 
around the polymer backbone is possible. The rescaled 
bending rigidity is given by the ratio of the persistence 
length and the rod length 
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(18) 



where L — 2aM is the polymer length. 
The external force can be caused by an external homo- 
geneous field, electric or gravitational. 
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where qe is the electric charge of the bead (g is the va- 
lency) for an electric field or the mass for a gravitational 
field as applicable to sedimentation experiments. The 
external force can also act on only part of the polymer, 
as for example the force applied by a molecular motor. 
Mutual penetration of monomers is prevented by an 
excluded volume interaction modeled by a truncated 
Lennard- Jones interaction 
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(20) 



'Tjl < 2 with an energy param- 



valid for separation jf^ 
eter a; = 3. 

Lateral periodic boundary conditions are implemented 
with the Lekner-Sperb summation scheme |70j for hy- 
drodynamic and electrostatic interactions. For suffi- 
cient numerical accuracy we choose time steps in the 



range jlo — 10~'^-10~®. Output values are calcu- 
lated every 10^-10^ steps, total simulation times are 
of the order of 10^-10^ steps, giving errorbars typically 
smaller than the symbol size. 



IV. HYDRODYNAMICS OF RIGID RODS 

We have seen that one of the principal feature of small 
Reynolds numbers hydrodynamics is that the fluid ve- 
locity field exhibits a linear dependence on exerted 
forces on immersed solid objects. In this paragraph, 
we explore this feature in some detail and demonstrate 
that propulsion at low Reynolds numbers is intimately 
related to shape design. The torque, assumed to act 
at the origin r == 0, and force exerted on a body are 
calculated by integrating the stress tensor a, Eq. ^, 
over the particle surface as [T] 



F = (t (T-dS 
N = * r X cr • dS 



(21) 
(22) 



It is plausible that the body shape - and thus any sym- 
metry properties of the solid - influence the flow at 
short distances, r > a. However for large distances 
r 3> a, we remind the reader that the flow field is 
invariably given by the Stokeslet, Eq. (Isll-Q. For a 
rigid object, the velocity v at any point r belonging to 
the solid is obtained by decomposing its motion into 
a global translation of velocity v'* and a rotation with 
angular velocity a;* related by v = v" + a;* x r. Using 
the no-slip condition at the body surface, one can then 
obtain the linear relationship between the forces F and 
torques N exerted on the rigid object and translational 
as well as rotational velocities v* and uj* by inverting 
equations (l2ll)-p2l) 



fJ-rt 



Mr 



(23) 



where the tensors /Xj and /x^ are the translation and ro- 
tation mobility tensors. Hence the off-diagonal tensors, 
/i.jr and /x^(, are the two coupling tensors and repre- 
sent cross-effects where rotation is induced by forces 
and translation induced by torques. Using the gener- 
alized reciprocal theorem [1], one can prove that the 
whole mobility tensor defined in Eq. ( 23 1 is symmet- 



ric. Moreover, if the solid object exhibits some type of 
symmetry, they will be reflected by symmetries of the 
mobility tensors /.i. 

In the following we are interested in bodies with three 
mutual orthogonal planes of symmetry (orthotropic 
bodies) such as ellipsoids, spheres, cylinders or cubes. 
In this case, we find /x^.^ = /x^^ = for every coordinate 
system. Hence translational and rotational motions are 
completely decoupled and only a torque can induce a 
rotation of the body. Moreover, tensors /x^ and /x^ are 
diagonal when expressed in the principal axes coordi- 
nate frame [1 . Hence for a cylinder, we have only two 



eigenvalues and the translation mobility tensor simpli- 
fies into 
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Mil 
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(24) 



where P is the rotation matrix into the principal axes 
coordinate frame. The two mobility constants /i|| and 
/U_L correspond to motions parallel and perpendicular 
to the cylinder, respectively. The exact values of these 
mobilities are only available for ellipsoids from which 
the values for a finite cylinder of length L have been 
approximately deduced. The literature gives the fol- 
lowing approximate values [T] 
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The important point is that for sufficient long cylinders 
these values can be legitimately approximated by 



Mil - 2mj 



(27) 



Suppose that a cylinder is falling under its gravitation- 
nal weight in a quiescent fluid at low Reynolds number. 
The angle /3 between its symmetry axis and the cylinder 
velocity v* is then given by (see Figure [2k) 
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where m is the cylinder mass, g the gravitationnal ac- 
celeration and a the angle between the symmetry axis 
and the horizontal (a' = 7r/2 — a). The deviation angle 
with respect to the vertical is given by [71 



a'-/3=--a-/3 
'2 



arctan 



tana' 
2 + tan^ a' 



(29) 



We thus recover the result that for a a cylinder parallel 
(a — 7r/2) or perpendicular (a = ) to gravity, it falls 
vertically (/3 — 0). The maximum fall angle is obtained 
for tan a' ~ \/2 and is around a' — /3 ~ 20° . 
Now we assume that the object is a circular helix with 
its axis pointing along z. We choose it arbitrary to 
be left-handed (see Figure [2J3). Planes perpendicular 
to the centerline are no more symmetric planes and 
the coupling mobility ^t^^ is non-zero. The rotation 
of the helix at angular velocity u) — ujz (in Fig. ^p 
we arbitraly choose w > 0) can now induce a thrust 
Ft = fJ-tr^^- A simple local argument of the generation 
of the thrust is given by H. Berg ^j and schematically 
shown in Fig. [2b. We approximate the helix filament 
by a chain of small cylinders, which make an angle S = 
2-KR/h. with the helix axis, where R is the helix radius 
and A the pitch. Because a cylinder has two different 
mobilites, Eq. (24), the viscous drag is twice as big 



when it moves sideways than when it moves end-on. 
Hence the components of the drag normal to the helix 
axis contribute to the torque whereas components of 
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FIG. 2: a) Oblique fall of a cylinder under gravity, b) 
Schematic analysis of the viscous drag on two cylinders of 
a helical filament contributing to thrust. The helix is left- 
handed and oj — ojz. The velocity of each segment is decom- 
posed into normal and parallel velocities, v± and V| | . Hence 
normal and tangential friction forces, Fj_ and _F||, act in op- 
posite direction and their ratio is F±/F\\ ~ 2v±/v\\. The 
total force has thus a component parallel to z, Ft- 



the drag parallel to the axis are all of the same sign and 
contribute to a propulsive thrust F^ — Ftz with Ft > 
0. Note that with our sign convention in Fig. [2b, we find 
fitr > 0. This sign is related to the helix handedness: a 
right-handed helix would give fj,tr < 0. The velocity is 
V = a; X R and we find for the tangential and normal 
forces to the filament per unit length 
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where p,.. and /ij^ are the friction cefficients per unit 



length. The thrust is then 



Ft = (-Fj_ sin 6 — F^^ cos S) z 

= -sin2(5 ( \ Ru: 

2 VM-L Mil. 



(32) 



Hence the thrust is maximized for d — 7r/4. By approx- 
imating Mil — 2/ij^, we finally find the coupling mobility 



Mr* 
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(33) 



Similarly, the rotation mobility associated with the re- 
sisting torque is given by 



Mr 



4/ij 



i?2(3 + cos2<5) 



(34) 



These approximate results allow a simple understand- 
ing of the propulsion mechanism. They are found us- 
ing the so-called resistive force theory developed by 
Gray-Hancock [73j and Lighthill ^EiW\- The roughest 
approximation is the neglect of hydrodynamic inter- 
actions between different segments which is of course 
not justified since these interactions are long-ranged 
in 1/r. Lighthill improves this approach by showing 
that these hydrodynamic interactions can be asymp- 
totically calculated in the limit of slender body the- 
ory [711 [7S], the small parameter being a/q where g <C A 
is the unknown typical length of small cylinders and a 
their thickness. Hydrodynamics are then taken into 
account by considering a helical line of Stokeslets for 
the far region (r > q) plus Doublets to fulfill the no- 
slip condition for the close region (r < q) where the 
helix is approximated by a straight filament. This al- 
lows a determination of q — 0.09A and local resistance 
coefficients per unit length p,7^^ = 2TTr]/\n{2q/a) and 

p.'J^ = 8TTi]/[2ln{2q/a) + 1]. Other works tackle the 
question of end effects and the influence of the cell body 
on the flow field, by considering Stokeslets, Doublets 
and a Rotlets which represents a point torque [76l [77] . 
Of course these studies suppose infinitely stiff helices. 



V. HYDRODYNAMICS OF SOFT RODS 

In the following, we consider elastic rods out of equi- 
librium under external forces. Contrary to rigid rods, 
hydrodynamic effects lead to a bending of elastics rods 
and therefore a reduction in symmetry [751 IZSj • It re- 
sults in an hydrodynamic translational-rotational cou- 
pling since now /i^j ^ 0, which can be samll but 
coupled to large forces and torques can lead to very 
different qualitative behaviour compared to infinitely 
rigid rods. Indeed, sedimenting rods orient perpendic- 
ularly to the direction of motion due to a hydrodynamic 
torque. Likewise, the bending of rotating rods can in- 
duce propulsion in the same way as helices, contrary 
to rotating stiff cylinders for which the frictional force 
is always perpendicular to the rotation axis. Hence a 
bent rod can play the role of propeller. 



A. Orientation of elastic rods under external field 

We consider an elastic rod submitted to an external 
force F'^^^ = qeE which can be the gravitational field in 
sedimentation experiments or an electric field (E) act- 
ing on a charged rod in birefringence experiments |15j . 
As briefly described in the introduction, the orientation 
of stiff charged rods such as tobacco-mosaic viruses in 
electric fields is mostly determined by polarizability ef- 
fects. The anisotropic electric polarizability favors in 
general an orientation with the direction of the largest 
polarizability parallel to the external field E. The 
largest contribution to the polarizability comes from 
the easily deformable counterion cloud close to each 
charged rod, which is maximal along the rod axis [5D1 - 
155] . As a result, at not too low fields, charged rod- like 



particles orient parallel to E. However, when hydrody- 
namic effects are dominant and overwhelm this induced 
dipole mechanism, rod like particles can orient perpen- 
dicular to the field E. In the following, we focus on this 
hydrodynamic orientation. 
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FIG. 3: Stationary deformation and orientation of an 
isotropic elastic rod at zero temperature, a) Snapshots for 
M = 19 and £*=0.03, 0.3 and 3 (from top to bottom) of a 
rod moving downwards, b) Rescaled mean curvature K/M 
(open symbols) and alignment parameter (x + 0.5)/M^ as 
a function of the rescaled rigidity e* confirming scaling pre- 
dictions Eq. (371 (broken lines with slopes -1 and -2). 



The mean chain bending is measured hy K = |foi — 
fAr_iAr|/2 = sin 6* where 9 is the bending angle of the 
terminal chain segment (see Fig. pa) and f are the nor- 
malized bond vectors. The overall rod orientation can 
be measured by the orientation parameter ip which is 
defined as 
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3 sin a ~ 1 



(35) 



where Re = r^v-i — Tq is the end-to-end vector of the 
chain, and the orientation angle a is the same as in 
Fig. [2^. The electric birefringence signal, on the other 
hand, is proportional to the alignment parameter |78j 



1 

M 



M-l 

E 

i=0 



E 



'u+i 



E 



(36) 



which is a measure of the average orientation of indi- 
vidual bonds and only equal to ip for straight rods. The 
hydrodynamic orientation mechanism, which of course 
only works at low Reynolds number, can be explained 



in simple terms as follows: the external force drives all 
monomers in the same way; the effective force, which 
results from the product of the hydrodynamic interac- 
tion tensor times direct forces, cf. Eq.(13l, is larger in 



the middle than at the two ends, because the middle 
receives hydrodynamic thrust from both sides. Since 
in stationary motion, all monomer velocities have to be 
the same, this imbalance in driving thrust is balanced 
by an elastic deformation, which in turn leads to an 
orientation of the rod. 

On a scaling level, the elastic torque of a rod of length 
L = 2aM bent by an angle 6, as shown in Fig. [3^, 
is Ng ~ e9/L and is balanced by an hydrodynamic 
torque, Nh, which arises from the above-mentioned in- 
homogeneity of hydrodynamic thrust. N^ is thus pro- 
portional to the direct force per bead, the number of 
monomers M and the rod length : Nh — L'^qeE/a. 
Equating Ng ~ Nh yields a stationary bending angle 
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(37) 



where e* = e/{qeEL'^) is the proper scaling variable in 
the continuum limit (a — >■ 0). In Fig. pp, we show the 
numerically determined rescaled chain bending K /M 
(open symbols) and the bond alignment (x + 1/2) /M"^ 
(dark symbols) as a function of the inverse driving force 
e* at zero temperature (the Langevin forces are switched 
off). Data are well described by the laws 6* ~ isT = 1.1 x 
lQ-^L/{ae*) and 6*2 ~ x + 1/2 = 1.1 x IQ-^L"^ /{ae*f 
(broken lines), for not too strong bending, in agreement 
with our scaling predictions. 

Rod bending reduces the symmetry and hydrodynam- 
ically couples translational and rotational degrees of 
freedom; due to a shift between the center of mass and 
hydrodynamic stress the bent rod is oriented perpen- 
dicularly to the direction of motion with the opening 
pointing backwards ^. For finite temperature the ori- 
entation is subject to thermal fluctuations and thus not 
complete. If we consider an assembly of elastic rods 
in solution as in sedimentation or birefringence exper- 
iments, we can deduce the average bending and ori- 
entation using Boltzmann statistics. Within linear re- 
sponse, the orienting torque N^ is proportional to the 
orientation angle a, the bending angle and the driving 
torque Nh'. N^ -~^ aONh- 

For low tem-peratures, the average orientational energy 
equals the thermal energy, aNa ~ ksT, leading to a 
mean square orientational fluctuation 
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where E is defined in Eq. ( 19 1. It is checked for a fixed 



rescaled chain stiffness e*L/a — (.^/{LE) in Fig. Wk 
with a numerically determined prefactor. The orienta- 
tional fluctuations at linear response are governed by a 
quadratic form proportional to a coupling constant J 
with expectation value 
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For low temperatures (high J-values) one obt ains 



1/J which together with Eq. (38) 



fixes the coupling constant J. 
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FIG. 4: a) Rescaled orientational parameter {ip + 0.5)M* 
as a function of the rescaled temperature kBT/{qeEa) = 
1/E for fixed rigidity e*L/a = 100 in agreement with the 
low-temperature scaling prediction Eq.pSl (broken line), 
b) Same data plotted as — ?/)/M* compared with the high- 
temperature scaling prediction Eq.(40l. 



The high-temperature regime is experimentally more 
relevant and corresponds to an almost isotropic rod ori- 
entation distribution. In this limit, one obtains from 
Eq. S (sin^a) ~ 1/3 - 2J/45 -|- ©(J^), which to- 



gether with the definition of ip, Eq. (35) yields the final 
result 



ih-^— = -1.1 X iO"^£;2— 

^ 15 4 



(40) 



-7r/2 



Equation ( 40 ) including the numerical prefactor is con- 
firmed in Fig. HJ3 (broken line) . As might be expected 
for the stationary behavior, the result is independent 
of the solvent viscosity. The dependence on the rod 
length -0 ~ L^ is stronger than for a mechanism based 
on electric counterion polarizability, where one typi- 
cally assumes an L-dependence like ij] ^ L^E^ or even 
smaller |80) . This means that hydrodynamic orienta- 
tion is dominant for long rods, in agreement with exper- 
imental observations showing increased anomalies for 
long polymers [TD]. Finally using Eq. (37) and (40 1, 
the ratio between chain orientation and chain bending 
becomes |'0|/^^ — ^p/L, which suggest that for stiff 
polymers £p > L one can observe weakly bent but still 
strongly oriented chains. 
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Transient birefringence and dichroism experiments [84] 
deal with the time dependence of the orientational re- 
sponse x{t) s-iid iplt) after the electric driving field, 
^ext _ gg^ jg suddenly turned off (or on) . The orienta- 
tional diffusion time of a rod, td = TuqeEfio/a = 218, 
is given by the perturbative expression for the rota- 
tional diffusion constant [5S] in powers of a/L. At the 
simplest level, the decay of the orientation angle a in 
an applied electric field follows a drift orientational dif- 
fusion and is governed by the equation 



da 
~dt 
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(41) 



where Na is the orientational torque. The hydrody- 
namic orientation time is therefore 
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hence within logarithmic corrections td/tyi^ — 1.1 x 
l{)-^L^{qeEf/{a^ekBT). The strong length depen- 
dence makes the orientational process much faster for 
long rods. For charged rod-like objects, the hydrody- 
namic orientation mechanism competes with the coun- 
terion polarization mechanism which favors parallel ori- 
entation. The orienting polarization torque is N^ ^ 
L^EQaE^ where £o is the dielectric constant (neglecting 
the dispute about the L dependence of the polarizabil- 
ity [10]). The characteristic orientation time due to 
electric polarization, Tp^ follows from a similar equa- 
tion as above, Tp ^ {kBT/D){a/N^) ~ l/(eoMoa-£'^): 
and is independent of the length. Hence, for long poly- 
mers one has Tp, > Tp > Thd, and the hydrodynamically 
induced anomalous birefringence is the fastest process. 
We illustrate the proposed mechanism with fd-viruses 
which have a length of L ~ 880 nm and a diameter of 
2a ~ 9 nm. The total net charge is roughly 500 e~ so 
that the valency per length 2a is about 9 = 5 [3]. For 
a typical electric field E = 10^ V/m the rescaled field 
strength is -E ~ 0.1. Equation (40 1 yields an orienta- 



tion of the order of unity when L/£p > 10 (a/L) ~ 
0.1, a realistic number for a virus material. In the 
case of sedimentation, the reduced driving field act- 
ing on the monomers \s E ~ Aiia^ pgG / [ikpT) where 
G = 9.81m/s^ is the gravitational acceleration, g is 
the g-factor of a centrifuge, and p ~ lO'^ kg/m'^ is the 
density difference between the sedimenting particle and 
the solvent. For a particle radius a ~ 10~^ m, we ob- 
tain E ~ 10~^5. To obtain the same effect as in the 
above charged-rod example, one would need a (7-factor 
of £( ~ 10^ which is large but reachable in an ultracen- 
trifuge. For larger colloidal rod-like particles, sedimen- 
tation will lead to hydrodynamic orientation even with- 
out centrifuging: large particles are predicted to sedi- 
ment perpendicularly to the direction of motion [13] . 
All our results are valid only for low Reynolds num- 
bers: for an object of length L moving at velocity m, 
the Reynolds number is Re ^ Lpu/rj. For a cylin- 
der the velocity scales as u ~ F/rjL where the total 
force is F ^ qeEL. The condition Re < 1 leads to 
{L/a){aqeE/kBT) ^ NE < Tf/{pkBT) w 10" (where 



we used the density and viscosity of water) which is ver- 
ified in our simulations. Hence the hydrodynamic orien- 
tation discussed here is compatible with low Reynolds 
number hydrodynamics. 



B. Propulsion with a rotating elastic rod 

We have seen in the previous section that the coupling 
between elasticity and hydrodynamics leads to surpris- 
ing effects: a straight elastic rod that is sedimenting 
in a quiescent fluid and thus subject to a homogeneous 
external force distribution becomes bent and oriented 
perpendicular to the force. Here we investigate a dif- 
ferent scenario: if one rotates an elastic rod by ap- 
plication of a torque at one of its ends, it will bend 
due to frictional forces. In the absence of bending, i.e 
for an infinitely stiff rod, a simple symmetry argument 
shows that the net thrust when averaged over a full 
turn vanishes. Clearly, a rod fixed at one end which 
rotates making a tilt angle with the rotation axis in 
the direction z, will induce no thrust along z. The 
cylinder moves sideways and the frictional force, given 
by F = p^_ijj X r, is always perpendicular to a; = wz. 
However, if the rotating rod is elastic, the friction forces 
will bend the rod and the question of propulsion thrust 
can be raised again. A few theoretical works have fo- 
cused on the coupling between hydrodynamics in a vis- 
cous medium and elasticity of soft polymers [86H88] . It 
was shown that finite stiffness of beating straight fila- 
ments breaks the time-reversal symmetry and enables 
propulsion. Self-propulsion using alternative mecha- 
nisms such as surface distortions [53] or chemical re- 
actions [Sn] has been actively proposed for designing 
macromolecular machines working in a viscous media. 
In this section, we present simulation results which in- 
dicate propulsion for rotating straight filaments [37] . 
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FIG. 5: Sketch in 3D (a) and projected views (b) of the 
rotating filament before (A, N = 600) and after (B, A*' = 
800) the bifurcation (M = 30 beads, i^jL = 10^). 



To mimic a rotary motor, we apply an external force 
on monomer 2 of a linear array of monomers, which is 
related to the applied torque N — Ni by F™' = N x 
ri2/r^2- The geometry and sign convention are shown 
in Fig. [5| Two different force ensembles are investi- 
gated: the stalled case, where the first two monomers 
forming the polymer base are held fixed in space by 
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applying virtual forces which exactly cancel all other 
elastic and hydrodynamic forces. In the moving case, 
we let the polymer move along the z direction and we 
apply virtual forces only laterally such that the base 
moves along a vertical straight track. A finite tilt angle 
at the base is imposed by a spontaneous curvature term 
in the elastic energy (e/2a)[l — cos(6'i — 91)]; here we 
show data for spontaneous curvature 91 = 45° . 
When a torque is applied to the filament base, it ro- 
tates and after a few turns exhibits a stationary shape. 
Due to friction the filament bends into a curved struc- 
ture. It should be noticed that, contrary to Section [V A[ 
the applied force is not homogeneous along the poly- 
mer. However the angular velocity of the whole poly- 
mer is the same for every monomer. Hence, the velocity 
V = a; X r, and therefore the frictional force, is larger 
at the end of the polymer, causing a stronger bending 
opposite to the direction of rotation. The bending due 
to hydrodynamic interactions (which are larger in the 
middle of the rod than at the ends), similar to effects 



studied in Section VA is thus a secondary and negli- 
gible effect. 

The stationary angular velocity uj = luA/iIq (in de- 
grees), as a function of the applied torque N — N/ksT^ 
is plotted in Fig. [6k. It shows a non-linear increase, 
which is associated to the finite bending rigidity. In- 
deed, for an infinitely stiff polymer, this relation would 
be strictly linear due to Stokes flow properties. For 
a critical torque, iVc, a shape bifurcation occurs and 
the angular frequency jumps dramatically. Figure [5] 
shows the stationary shape of the polymer before (A) 
and after the bifurcation (B). The high-torque state is 
characterized by a smaller distance r_L from the rotation 
axis. This explains the observed jump in angular veloc- 
ity since the applied torque, i.e. the frictionnal forces, 
remain constant at the bifurcation. Since the velocity 
of a bead z is v; = a; x r^, we find uj = Vi/r±^i which 
jumps dramatically since rj_i decreases. Pronounced 
hysteresis, which becomes slightly weaker with decreas- 
ing torque sweep rate, is observed when N is varied 
across the critical region, as illustrated in Fig. |6]d. This 
hysteresis cycle indicates that this shape bifuraction is 
subcritical. 

What is the control parameter of this bifurcation? Fig- 
ure |6j: shows the dependence of the critical torque on 
the persistence length of the rod. Since the torque 
sweep rate is finite, increasing (o) and decreasing torque 
(D) give rise to slightly different bifurcation values. For 
very small persistence lengths, conformational fluctu- 
ations suppress the transition. The shape bifurcation 
can be understood by simply balancing elastic and driv- 
ing torques: the bending torque due to the filament 
deformation, projected along the vertical axis, reads 
~ £sm9i/R where R is the bending radius of the fila- 
ment, which has to be balanced by the external torque 
N . The onset of the transition is fixed when the bend- 
ing radius reaches the length of the filament, R ^ L. 
This yields a critical torque 
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FIG. 6: a) Angular velocity as function of external torque, 
N, for various (.^/L: 3333 (o), 6667 (D), 1.33 x 10^ (o) and 
2x 10*(o) [inset: 333 (o), 1000 (D) and 2500 (o)]. b) Hystere- 
sis cycle for l-p/ L — 1667 for increasing (o) and decreasing 
torque (D) at constant rate SN/{St/A) = 1.25 x lO"**. c) 
Persistence length vs. critical torque for increasing (o) and 
decreasing torque (D) both obeying a linear law according 
to Eq. ((43| (AH data for stalled case, 14 = 0). 



which is shown in Fig. [6]; as a solid line and agrees 



very well with the numerical data. It transpires that 
the bifurcation is purely elastic in origin, since nei- 
ther hydrodynamical parameters nor temperature ap- 
pear in Eq. ( 43 ) . Using the rotational mobility ^,.ot ~ 
l/(ry_L'^), the lower critical angular frequency reads 
'^c ^ A'rotA^c ~ esm9i/{r]L'^). This threshold turns 
out to be much lower than for the continuous twirling- 
whirling transition of a rotating rod with torsional elas- 
ticity, wj^ ^ £ / [rjo? L"^) ^ which was obtained using lin- 
ear analysis |88j . Moreover, the twirling-whirling tran- 
sition is related to a symmetry breaking, named as su- 
percritical bifurcation. 
At this point, one may wonder what the role of hy- 
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drodynamic interactions is in this transition. This role 
is twofold: first, it induces a propulsive thrust, since 
we have seen in Section |IV[ that hydrodynamic inter- 
actions are essential in the production of thrust by a 
rotating helix. Obviously, our bent filament for large 
torques does not have the geometry of a pure helix but 
it wraps around the rotation axis. The second aspect is 
more subtle. To understand the role of hydrodynamic 
interactions for the shape bifurcation, we implemented 
two additional approximations: 

• the free draining approximation (FDA), where 
hydrodynamics are ignored altogether. The mo- 
bility matrix is thus isotropic: /i^ = /iq 1 ^ij] 

• the resistive force theory within the slender body 
approximation (SBA) which is defined in Sec- 
tion |IV| |18| . The mobility matrix is given by 
Eq. ( |24p . The choice of the values of the two local 
mobilities in the simulations is of course difficult 
and heuristic, since at our coarse-grained scale the 
monomers are spheres which sets fj,± — fJ-u — fJ-o- 
Values for long cylinders or long ellipsoids are 
given in literature and we arbitrarily chose fi± — 
3no{2\nM + l)/4 and ^^ = 3^io{2\nM - l)/2. 
Note that these values should in principle be ad- 
justed to the shape of the polymer. 

Figure [7]^shows the angular velocity versus the applied 
torque, N, for the three approximations. In accord with 
our elastic scaling model, the transition occurs for all 
approximations even in the absence of hydrodynamics 
(FDA: o). However, the critical torque is overestimated 
for both FD and SB approximations by roughly a factor 
of 5. 
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FIG. 7: Angular velocity vs. applied torque for i!p/L — 1667 
using full hydrodynamic interactions ( x ) , slender-body (D) 
and free-draining approximations (o). The transition occurs 
in all cases but hydrodynamics decrease significantly the 
critical torque (Stalled case, 14 = 0). 



Is it possible to define the mobilities of the whole ob- 
ject? The relation between angular velocity and torque 
is no more linear and the meaning of the matrix in 
Eq. ( 23 ) changes. We define four mobilities of the whole 
object following Eq. (23). The dimensionality of the 



mobility matrix is now 2x2 since we consider only the 
z-component of torques, forces and velocity of the whole 
rod: V = V^z and F = F^xtz the corresponding exter- 
nal force applied at the propeller base. Nevertheless, we 



shall keep in mind that the matrix notation is just used 
for an easier comparison to rigid propellers, because 
these four mobilities (/Zi-r, /Zj-t, Mtr and /Xtt) depend a 
priori on N and F°'^^. Moreover, for a flexible rod the 
symmetry of the mobility tensor is not satisfled. The 
propulsion velocity along the rotation axis is plotted in 
Fig. 8a as a function of the control parameter NL/ip 
for dmerent persistence lengths for the z-moving case, 
i.e with i^ext = 0. At the transition {NL/lp ~ 0.7), a 
jump in the propulsion velocity is observed and Vz is 
almost linear for N > N^, revealing that the shape re- 
mains almost fixed in this range of torque values. Inset 
of Fig. [8^ shows the ratio of 14 before and after the 
transition against the same ratio for angular velocities. 
The variation is roughly linear, meaning that the jump 
in propulsion is mostly due to the increase in angular 
velocity. A refined analysis would be necessary to de- 
tect if the rod takes a specific shape which favors the 
forward velocity. For instance, in the case of a per- 
fect helix, we have seen that the thrust is maximised 
for a pitch angle equal to a = 45°. We have checked 
that all the previous results for uj{N) are very similar 
in the z-moving case. The only difference appears at 
the bifurcation, the critical torque being a little smaller 
without any external load force. 

Up to now, our model neglected the presence of a base 
for the rotor. In reality, the filament is attached at its 
end to a body, the head of the whole propeller. We 
assume that the relative rotation of both components 
is free. The net force and torque on the complete pro- 
peller (body -|- filament) swimming at constant velocity 
are zero |16| since the motion is overdamped preventing 
any acceleration. Hence, for a spherical body of radius 
R and angular velocity fl, the loading force is propor- 
tional to the body velocity whereas the loading torque 
is proportionnal to J7, which reads 



-fcxt 

N 



-67rr]R 14 
-87rryi?3 fl 



(44) 
(45) 



To describe completely the whole propeller with a per- 
fect load matching, it is thus necessary to find the char- 
acteristic response of the filament to an external load. 
To test numerically our propeller under load, we ap- 
plied an external force, -Foxt: which we define to be 
positive when it pushes against its natural swimming 
direction. Figure [Sp shows the variation of Vz with 
Fcxt for £p/L = 1667 at two different torques, just be- 
low and above the bifurcation. The laws are almost 
linear in both cases, meaning that ^tt is almost inde- 
pendent of Fext- Hence for given torque N and body 
radius R, load matching fixes i^ext and thus T4 by tak- 
ing the intersection between the above law and Eq. ( 44 1 . 



The angular velocity of the body is then found using 
Eq. ( 45 1 . The stalled case corresponds to an infinite 
spherical body (R — >■ oo) whereas the free z-moving 
case imposes R = 0, i.e. no body. Of course, our study 
is in principle not valid for i? 3> 1 because the presence 
of a hard wall will perturb the fiow at the filament base. 
The efficiency of the power converter can be defined as 
the ratio of the propulsive power output and the rotary 
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FIG. 8: a) Propulsion velocity parallel to the rotation axis 
vs. NL/£p in the z-moving case for (.-p/L = 2500 (x), 5000 
(D) and 2 x 10* (o). The inset shows the variations of the 
ratio of velocities after and before the transtion vs. the 
ratio of angular velocities (o, l-p/L = 1667). b) Vz vs. the 
external force _Fcxt applied on monomer 2 {(-p/ L = 1667) 
before the bifurcation for iV = 0.787Vc (D) and after iV = 
0.84iVc (o). c) Efficiency, 77, vs. Fcxt (same parameters as 
b). The solid line is a parabolic fit. 



power input, 



-Fcxt Vz 

Nuj 



By inserting Eq. (23) in Eq. (46), we obtain 



ri{F,.t) 



Mtr^^cxt 



Mtt-Fexl 



AirrA^^ + MrtA^i^c, 



(46) 



(47) 



We have checked that /irt is negligibly small. Hence, the 
efficiency becomes parabolic as a function of the exter- 
nal force as shown in Fig. |8b. The highest efficiency 



is obtained for Fcxt = ~^trN/{2fitt) = -Fstaii/2 and is 
only of the order of 1% after the transition and 3 times 
smaller before. With this load, we find a body size 
i? ~ 5 a and an unique \n\ = l{a/R)^N ~ 8 < w ~ 250 

for A^ = 0.84A'c after the transition. 
In conclusion, we show that a simple straight elastic fil- 
ament can induce propulsion when rotated at one end. 
The interplay between elastic deformations and hydro- 
dynamic interactions results in a substantial directed 
thrust due to a subcritical dynamic bifurcation for a 
control parameter N/Nc larger than 1. The mobility 
/itr changes its sign when the torque is reversed, which 
implies a forward thrust whatever the sense of rotation. 
This work provides a clue for the synthetic manufac- 
ture of biomimetic micro-propeller, by using simple 
semi-ficxible polymers instead of rigid helices. This 
bifurcation could be experimentally investigated using 
a macroscopic scale model similar to the one devel- 
oped by Powers group 9.1]. With their experimental 
setup, the bifurcation would occur for torques on the 
order of 0.05 N.m and angular velocities larger than 
0.01 Hz, which are accessible. Another type of exper- 
iments has been developed by Bibette's group, where 
they make microscopic artificial swimmers with chains 
of connected ferrofluid droplets swimming by applying 
a transverse oscillatory magnetic field |SJ|. The result- 
ing forward thrust is due to a beating of this "elastic" 
chain. Even if these experimental results are associ- 
ated to transverse beating, one can imagine quite soon 
a rotating chain experiment. 



C. Grafted elastic polymers in shear flow 

In previous theoretical studies, the effects of shear 
flows on surface-anchored polymers have been primarily 
studied analytically and numerically for dense brushes 
and flexible polymers [55ti^ . In this section, we study 
the shear flow response of a layer of semiflexible poly- 
mers grafted to a planar surface (a so-called polymer 
brush) as a function of polymer stiffness, grafting den- 
sity and shear rate f96|. A finite orientational stiffness 
is introduced at the grafting points, tending to ver- 
tically orient the grafted chains at the basal region. 
The time evolution of the system is described by the 
Langevin equation ( 13 1 where a linear external shear 
flow v(r) 



7ZX along the x direction with a shear rate 
7 is included. We explicitly incorporate full hydrody- 
namic interactions between monomers in the presence 
of a no-slip surface through a non-diagonal mobility 
tensor fiij on the Rotne-Prager level [3j following Sec- 
tion |irc] As before, we rescale all lengths and energies 
and express the rescaled shear rate by 7 = a^jij/ksT. 
In order to mimic a semi-inflnite surface covered with 
anchored chains at a finite density, we consider 4 poly- 
mers with monomer number A'^ = 10 grafted regularly 
within a square unit cell of lateral size D and treat 
hydrodynamic interactions via lateral periodic bound- 
ary conditions imposed on the flow fleld. The grafting 
density is defined as p = Aa^ /D^. 
The effects of a shear flow is demonstrated in typical 
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FIG. 9: Snapshots of polymer configurations for grafting 
density p = 0.005 and semiflexible polymers (top row, per- 
sistence length ip/L = l)and rather stiff polymers (bottom 
row, persistence length Ip/L = 5). To the left no shear flow 
is acting, in the middle column the shear rate is 7 = 0.003 
(corresponding to a shear rate in unrescaled units of 7 ~ lO'^ 
s~^ for sphere radius a — Inm and water as a solvent), and 
to the right the shear rate is 7 = 0.03 (corresponding to a 
shear rate in unrescaled units oi j ^ 10* s~^ for a = 1 nm). 



simulation snapshots shown in Fig. |9] The upper row is 
for semiflexible polymers, while the bottom row shows 
results for rather stiff polymers. As one goes from the 
left to the right, the applied shear rate is increasing 
steadily, and one observes a clear change of the polymer 
conformations. Upon shearing the brush, the polymers 
tend to be bent to reduce hydrodynamic friction, and in 
turn exert a resistive response to the flow due to a flnite 
bending stiffness. Therefore, the resulting solvent flow 
velocity is expected to be reduced compared to the case 
of a bare surface in the absence of grafted chains. As 
the shear rate increases, the bending of the polymers 
becomes stronger and they move closer to the surface. 
As a result, the flow can approach the no-slip surface 
more closely. 

Figure [TO] shows the stationary solvent velocity pro- 
files v{z) = Av/{12tt'^L) measured from the simula- 
tion, and clearly exhibits a shear-dependent lift-up of 
the shear plane. By linear extrapolation of the flow 
profile from large distances to the surface, the shear- 
dependent stagnation thickness S (corresponding to the 
lift-up of the shear plane) in units of a is quantitatively 
determined in the figure (dashed lines). It is numer- 
ically difficult to reach the small shear rates typically 
used in experiments. We instead perform simulations 
for elevated shear rates at which numerical errors are 
small and extrapolate down to small shear rates. Let us 
discuss the experimental relevance of parameters used 
in the simulation: we rescaled every parameter using 
monomer radius a and thermal energy fc^T. There- 
fore, choosing a proper value of a, we are able to cover 
a wide range of experimental situations. For example, 
when a = 1 nm is assumed, as applicable to DNA, 
it gives a contour length L — 20 nm and shear rate 
7 fa 10* s~^ for 7 = 0.03 (which is a quite big number 
in a conventional set-up). However, assuming a = 100 
nm, one gets L = 2/xm and 7 « 10^ s~^ for 7 = 0.03, 
which can be easily reached in the laboratory. 
In order to compare with the experimental results, 
we study the stagnation thickness for different applied 



V 




FIG. 10: Solvent flow profiles as a function of distance from 
the surface for grafting densities p — 0.003 (diamonds) and 
p = 0.007 (triangles). Parameters of the system are chosen 
as (.pjL = 5 and 7 — 0.03 (which corresponds to 7 ~ 10* s~^ 
for a = 1 nm). The length of polymers is roughly L — 20a. 
The effective shear-dependent stagnation length S is deter- 
mined by a linear extrapolation of the flow profile towards 
the grafting surface (dashed lines) . For comparison, the un- 
perturbed flow velocity in the absence of grafted chains is 
denoted by a solid line. 
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FIG. 11: Stagnation length 5 vs. shear rate 7 for two dif- 
ferent grafting densities: p = 0.005 (filled diamonds) and 
0.001 (open diamonds). The unitless shear rate 7 actually 
used in the simulation has been converted into an explicit 
shear rate 7 via 7 — kBTj/a'^ri and assuming a monomer 
radius of a = 1 nm, as applicable to DNA. The persistence 
of the filament is £p/L = 5. 



shear rates. As shown in Fig. [TT] the stagnation length 
6 decreases for increasing shear rate. While the chains 
are strongly bent close to the surface under strong 
shear, they tend to keep an upright conformation under 
weak shearing, leading to significantly enhanced hydro- 
dynamic friction. 

This is in qualitative agreement with recent laser- 
trap experiments on DNA chains that are anchored on 
spherical colloids [28] which results are reproduced in 
Fig. [12] In the experiments the brush-covered colloids, 
which have a diameter of roughly a micrometer, are 
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FIG. 12: Ratio of the efTective hydrodynamics radii th for 
DNA-grafted and blank colloids vs. flow velocity. Repro- 
duced from [28| . The symbols corresponds to three col- 
loids with different DNA lengths (o: 1000 bp; A: 4000 bp; 
D: 6000 bp and grafting density of 0.03 ± O.Ol/iTn^) in two 
NaCl electrolytes (open symbols: 150 mM; filled symbols: 
100 mM. The results show that the hydrodynamic radius 
decreases with increasing shear rate, indicative of confor- 
mational changes of the grafted polymers and an effective 
decrease of the brush layer. 



held fixed in tlie laser trap and subjected to laminar 
flows of the order of up to 1000 ^m/s, giving rise to 
shear rates of the order of 1000 s~^. The force acting 
on the colloid can be measured and thus the hydro- 
dynamic radius in inferred. We observe that the hy- 
drodynamic radius of DNA-grafted colloids decreases 
with the flow velocity, i.e. the shear rate in our model. 
Hence, the stagnation length decrease with increasing 
shear rate, which constitutes a non-linear friction re- 
sponse of the suspended colloidal sphere as a function 
of solvent velocity. This is in qualitative agreement 
with the simulations. 



teractions and elastic deformations of soft matter plays 
an important role. When a significant external force 
is applied to the monomers of a deformable polymer, 
a complex hydrodynamic interaction mediated by the 
surrounding water sets in and couples the velocities 
of different components of the system. This can lead 
to unexpected motions of the whole system: i) a 
sedimenting elastic rod deforms and, as a consequence, 
aligns perpendicularly to the direction of motion; ii) 
a rotating filament deforms and therefore gives rise 
to propulsion perpendicularly to the rotation axis, 
and Hi) the deformation of surface anchored polymers 
in shear flow leads to a shear-rate dependent shift 
of the stagnation-plane position. Since all materials 
have actually finite elastic moduli, their deformations 
in flow fields are unavoidable. These phenomena 
are therefore relevant to recent experiments probing 
the coupling of elasticity and hydrodynamics of 
nano-systems far from equilibrium. It is suggested 
that hydrodynamic interactions should also be taken 
into account in biological systems [87]. Indeed, as an 
example, one may consider the propulsion mechanism 
studied in Section [VB| and carry the lesson to bacterial 
propulsion. Of course, flagella are helical and their 
physics is therefore much more involved than a simple 
straight propeller. However, flagellar motors generate 
sufficient torques for the non-linear elastic phenomena 
discussed in this review to occur: they are powered by 
a proton-motive force which yields torques on the order 
of iV « 10^ keT 123 [SS]. Moreover, the flagellum 
length is i ~ 10 ^j,m, leading to £p/L ~ 10^ - 10* 
and NL/£p ~ 0.1 — 1. Hence, torque-induced shape 
transformations are most likely biologically relevant 
and might be directly observed with straight biopoly- 
mers attached to flagellar motors. Indeed, the role 
of flexibility in bacterial propulsion has been recently 
stressed since bending is crucial for the bundling of 
flagella [911[99]- 



VI. CONCLUSION 

Three different situations have been reviewed having in 
common that the coupling between hydrodynamic in- 
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